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ABSTRACT 

We show how hydrostatic equilibrium in galaxy clusters can be quantitatively probed combining X-ray, SZ, and gravitational-lensing 
data. Our previously published method for recovering three-dimensional cluster gas distributions avoids the assumption of hydrostatic 
equilibrium. Independent reconstructions of cumulative total-mass profiles can then be obtained from the gas distribution, assuming 
hydrostatic equilibrium, and from gravitational lensing, neglecting it. Hydrostatic equilibrium can then be quantified comparing the 
two. We describe this procedure in detail and show that it performs well on progressively realistic synthetic data. An application to a 
cluster merger demonstrates how hydrostatic equilibrium is violated and restored as the merger proceeds. 
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1. Introduction 

Numerous observations show that galaxy clusters frequently ex- 
hibit irregular shapes and violent dynamics. On the theoret- 
ical front, simulations indicate that cluster-sized dark matter 
halos are often well described as triaxial ellipsoids, but not 
as spheres ( |Jing & Suto||2QQ2] ). Nevertheless, clusters are of- 
ten interpreted as spherically- symmetric objects in hydrostatic 
equilibrium, which is a potential source of error. For example, 
[Hallman et al.| ( |2006 | show that intrinsic variations in clusters 
limit the accuracy of cluster gas mass estimates to about 10% 
when using such simple assumptions. 

Several authors tried to relax the restricting assumption of 
spherical symmetry a nd aimed at a joint analysis of different 
types of cluster data. |Zaroubi et al.| ([1998 ) suggested to base 
the reconstruction of axisymmetric, three-dimensional gravita- 
tional cluster potentials on the Fourier slice theorem, extrapolat- 
ing Fourier modes into the "cone of ignorance". They applied 
their technique to simulat e d data and s howed that it performs 
well HZaro ubi et "^1)120011). |Dore et al. (2001) followed a per- 
turbative approach, and Lee & Suto (2004 ) proposed to adapt 
parameters of triaxial halo models, all by combining different 
data sets such as X-ray, (thermal) Sunyaev-Zel'dovich (SZ) and 
gravitational-lensing maps. A similar method was applied to data 
by De Filipp is et al. (2005) . 

An alternative approach based on the iterative Richardson- 
Lucy deconvolution was suggested by |Reblinsky| ([2000) and 
[Reblinsky & Bartelmann (2001). It aims at the gravitational 
potential, assumes only axial symmetry of the main cluster 
body, avoids extrapolations in Fourier space, a nd can easily 
be extended to include additional data sets. In IPuchwein &I 



Bartelmann (2006) (hereafter Paper I) we developed the latter 



algorithm further. However, instead of aiming at the gravita- 
tional potential, which would require us to assume a relation 
between the gas distribution and the gravitational field, we pro- 
posed methods to reconstruct the three-dimensional cluster gas 
density and temperature distribution from X-ray and thermal SZ 
effect observations. These methods do not require any equilib- 
rium assumption other than local thermal equilibrium and again 
assume only axial symmetry with respect to an arbitrarily in- 



clined axis. Using synthetic observations of analytically mod- 
elled and numerically simulated galaxy clusters we showed that 
these reconstruction methods perform very well, even in the 
presence of observational noise, deviations from axial symme- 
try and cluster substructure. 

In this work we use this gas reconstruction algorithm to- 
gether with novel methods to reconstruct the three-dimensional 
gravitational potential from lensing data in order to probe hy- 
drostatic equilibrium in galaxy clusters and to quantify the ac- 
curacy of mass estimates based on the assumption of hydro- 
static equilibrium. We will also introduce methods to find three- 
dimensional reconstructions of the gravitational potential and of 
the mass profiles of relaxed galaxy clusters from X-ray and ther- 
mal SZ observations alone. All these methods are tested with 
synthetic observations of analytically modelled and numerically 
simulated galaxy clusters. Mass estimates based on the gas re- 
construction and the assumption of hydrostatic equilibrium are 
compared to lensing mass estimates and to the original analytic 
or simulated masses. 



2. Three-dimensional cluster reconstruction 
techniques 

2.1. Cluster gas reconstructions from X-ray and SZ data 

In Paper I a novel technique to reconstruct the intra-cluster 
medium in three dimensions by a combined analysis of X-ray 
and thermal Sunyaev-Zel'dovich effect observations was intro- 
duced (hereafter also called XSZ reconstructions). It assumes 
only axial symmetry of the cluster halo with respect to an ar- 
bitrarily inclined axis and does not require any equilibrium as- 
sumption other than local thermal equilibrium. The iterative 



method is based on Richardson-Lucy deconvolution ( [Lucy 1974 



,1994 ) and is a generalisation of employing it to reconstructions 
of three-dimensional axisymmetric quantities form their projec- 
tion along the line-of- sight ( Binney et al.| 1990). Due to the as- 
sumed axial symmetry the reconstructed gas densities and tem- 
peratures depend only on the distance R from the symmetry axis 
and the coordinate Z along the axis. 
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In this study we employ the same three-dimensional gas re- 
construction method as in Paper I except for one modification. 
We now use a more realistic model for the X-ray emission of the 
intra-cluster medium that also includes line emission. The syn- 
thetic X-ray observations of the reconstructed cluster halo that 
are performed during the iterative deprojection (see Paper I) are 
now calculated with the MEKAL emission model (see Liedahl" 
|etaL] ( [T995] ); |Kaastra & Mewe (1993)) and the WAR S model for 
galactic absorption ( [Morrison & McCammon|1983| ). More pre- 
cisely we use th e X-ray spectral fitting software package XSPEC 
( |Arnaud|1996[ ) to create a table of the cooling function with the 
models mentioned above, assuming a constant metallicity of 0.3 
times the Solar value and an equivalent hydrogen column den- 
sity of 5 X 10^^ atoms cm~^. We then use this table to produce 
the synthetic X-ray maps needed during the reconstruction and 
we also employed it to create the synthetic observations of an- 
alytically modelled and numerically simulated clusters that are 
discussed in sections [3] and |4] Note that it is not necessary to 
change the equations for the iterative corrections of the gas den- 
sity and temperature (Eq. (22) and (23) in Paper I), which were 
derived considering only thermal bremsstrahlung, because small 
errors introduced by using these equations are corrected in sub- 
sequent iteration steps. It also does not significantly affect the 
number of iterations needed to achieve a good reconstruction. 

2.2. Reconstructions of the gravitational potential and the 
total mass distribution from X-ray and SZ data 

We want to use these three-dimensional cluster gas reconstruc- 
tions to find the gravitational potentials and total mass distribu- 
tions of relaxed galaxy clusters and to probe the dynamical state 
of potentially unrelaxed clusters by comparing such gas recon- 
structions to an analysis of lensing data. 

To find the gravitational potential of a cluster from the distri- 
bution of the cluster gas we assume that the gas is in hydrostatic 
equilibrium. Then the gas density p, the gas pressure p and the 
gravitational potential (j) satisfy 



V0 : 



p ' 



(1) 



In principle this equation can be used to find the gravitational 
potential of relaxed clusters from three-dimensional reconstruc- 
tions of their intra-cluster medium. However due to deviations 
from hydrostatic equilibrium, and the presence of observational 
noise and cluster substructure violating axial symmetry, the curl 
of —Vp/p will not vanish exactly for the reconstructed gas dis- 
tributions. Thus one cannot obtain a unique solution for di- 
rectly from Eq. ([T]). 

To get a unique solution we first derive —Vp/p on the grid 
in R and Z space on which the gas reconstruction was calculated 
(see also Paper I). Then we aim to determine the potential for 
which V0 is closest to — V/?/p. We do that by finding the values 
of the potential at all grid points which minimise the deviation 



neighbours ij 



Pj-Pi 



-PiY 



(2) 



between these two vector fields. Here Pi^Pj are the gas pres- 
sures, pi^pj the gas densities and 0;, 0j the gravitational poten- 
tials at the R and Z coordinates of grid points / and j. The sum 
extends only over such pairs of grid points / and j that are near- 
est neighbours. Conjugate gradient minimisation starting with a 
guess = is used to find the solution for the 0^ . 



However, to reduce noise in the potential it turned out to be 
favourable to add a penalty function to ^ that requires the sec- 
ond derivatives of the potential to be small. We also multiply 
each term in the sum in Eq. ^ and the penalty function by a 
weight factor. So the function we end up minimising is. 



neighbours ij 



\{pj+pi)' 



- (^iR^ - 200' + (te> + te< - 200') , (3) 



where rij is the distance from the cluster centre to the midpoint 
of the line connecting grid points / and j. /7?>, /7?< and iZy.JZ^ 
are the indices of the neighbouring grid points of point / in the 
R and Z directions respectively. The weighting function w(r) 
is chosen equal to one in the central region of the cluster, for 
r < 0.3/, then it smoothly goes to zero, and vanishes for r > 0.4/, 
close to the perimeter of the box with side length / that is used for 
the gas reconstruction (see Paper I). This is necessary because 
there are significant artefacts in the gas reconstruction close to 
the perimeter (see also Paper I). When using the potential recon- 
struction algorithm proposed here they would have a non-local 
effect on the potential reconstruction and would thereby reduce 
its quality also near the cluster centre. The weight factor Wp = 3 
for the penalty function was chosen by trial and error and proved 
to be effective. 

One can then use the reconstructed three-dimensional gravi- 
tational potential to find the total mass distribution of the galaxy 
cluster. 

A simpler alternative way to get mass estimates from the gas 
reconstruction and the assumption of hydrostatic equilibrium, 
that does not require us to reconstruct the gravitational poten- 
tial, is to apply Gauss's law to the gravitational field and use Eq. 
([T]) to express the gravitational field as V/7/p. This allows us to 
define a cumulative mass M<^ xsz as a function of radius r from 
the cluster centre by 



1 f Vp 



(4) 



where G is Newton's constant and the integral extends over 
the surface of a sphere with radius r around the cluster centre. 
The numerical evaluation of the integral is done using 128 sam- 
pling points which are equally spaced in the angular coordinate 
= arctan(Z/7?). For each point the component of Vp/p per- 
pendicular to the surface is calculated from the gas reconstruc- 
tion and multiplied by the area of the corresponding ring. 

2.3. Reconstructions of the gravitational potential and the 
total mass distribution from lensing data 

Lensing obs ervations allow reconstructions of the lensing poten- 
]('2006 )), which is simply the suitably 



tial (see e.g. Cacciato et al. 



rescaled projection of the lens gravitational potential along the 
line-of-sight. Once the lensing potential is found, Richardson- 
Lucy deconvolution can be applied to deproject it in order to 
obtain the three-dimensional gravitational potential. Again axial 
symmetry with respect to an arbitrarily inclined axis needs to be 
assumed. 

We employ the deprojection algorithm discussed in sections 
2.1 and 2.2 of Paper I to obtain such three-dimensional recon- 
structions of the gravitational potential. 

In Paper I the optimal number of iterations for three- 
dimensional gas reconstructions was studied. Richardson-Lucy 
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deconvolution reproduces large scale structure quickly, while it 
converges slowly to small scale structure. It turned out that for 
gas reconstructions based on X-ray and SZ data it is best to use 
about five iterations. For a smaller number of iterations the clus- 
ter structure is not recovered sufficiently well, while for a larger 
number of iterations the reconstruction algorithm tries to repro- 
duce small-scale observational noise which can reduce the re- 
construction quality again. However as the lensing potential is a 
much smoother quantity than the X-ray surface brightness or the 
SZ temperature decrement it is favourable to use a larger number 
of iterations for deprojections of the lensing potential. 

However even when using a large number of iterations, prob- 
lems with the gravitational potential reconstruction arise for 
small inclination angles / between the line-of- sight and the sym- 
metry axis, because then the assumption of axial symmetry con- 
tains least information (see also Paper I) and a reconstruction 
that, compared to the original halo, is stretched along the sym- 
metry axis can still reproduce the lensing observations rather 
well. For a cluster with a roughly spherical gravitational poten- 
tial and for a small inclination angle one gets too large correction 
factors close to the symmetry axis during the first few iteration 
steps when starting from a flat guess and thus the reconstruction 
after a few iterations is overly extended along that axis. As the 
power to determine the halo elongation along the symmetry axis 
is limited for small inclination angles the reconstruction algo- 
rithm takes very long to recover from this. To avoid this problem, 
it is thus favourable to start with a guess that has already more 
or less the right shape. We can get such a guess by doing a grav- 
itational potential reconstruction from a flat guess with a small 
number of iterations and by then making the obtained potential 
spherically symmetric while preserving its profile. This spheri- 
cally symmetrised potential can then be used as a first guess for 
the actual reconstruction with a larger number of iterations. 10 
iterations were used to produce spherically symmetrised guesses 
for reconstructions from synthetic lensing data in sections [3] and 
m The actual reconstructions use 30 iterations and start either 
from such a spherically symmetrised guess or from a flat guess 
as specified there. 

The lensing three-dimensional gravitational potential recon- 
structions can then be used to find the total mass distribution and 
can be compared to reconstructions from X-ray and SZ data. In 
order to have a quantity that can be directly compared to M<^ xsz 
we define in analogy to Eq. ^ a lensing cumulative mass 

1 
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where is the three-dimensional gravitational potential obtained 
by deprojecting the lensing potential. The numerical evaluation 
of the integral is done in the same way as for M<^ xsz- 

2.4. Probing hydrostatic equilibrium 

Hydrostatic equilibrium in galaxy clusters can be probed by 
comparing cluster reconstructions based on X-ray and SZ data 
to lensing reconstructions. In principle this could be done by 
comparing the gravitational potential obtained by minimising 
Eq. ([3]) to the one found by deprojecting the lensing potential. 
However as the gravitational potential is not uniquely defined it 
is more favourable to compare the cumulative masses defined in 
Eqs. ^ and ([5]). If the cluster is exactly in hydrostatic equilib- 
rium, so that Eq. ([T]) is satisfied, the masses should be identical 
for all distances r from the cluster centre except for small de- 
viations caused by reconstruction errors. Otherwise differences 
between the masses directly reflect the differences between the 



gravitational field and V/7/p. In the next two sections we test 
this method to probe hydrostatic equilibrium in galaxy clusters 
by performing such a comparison using synthetic observation of 
analytically modelled and numerically simulated clusters. 

3. Probing hydrostatic equilibrium in analytically 
modelled clusters 

3.1. The analytic halo model 

We use an analytic halo model with a NEW total (gas+DM) den- 
sity profile to test the methods introduced in section [2] Thus the 
total matter density pm and the gravitational potential are given 
by 

c 

^ = ^lnm, (7) 



where is the NEW scaling radius and c = 4p(r^) fixes the nor- 
malisation of the density profile. Eor the cluster gas we assume 
in this toy model that the ratio / of V/?/p to — V0 is constant but 
can be different from 1 . So Eq. ([T]) generalises to 



/V0 = - 



(8) 



We further assume a poly tropic equation of state T oc pl'~i for 
the cluster gas, where T is the gas temperature and 7 the poly- 
tropic index. Then the gas density p and temperature T satisfy 



(i-r)<^> 



kT = f m, 



(9) 
(10) 



where k is Boltzmann's constant and m is the mean gas particle 
mass. In the following, we adopt 7= 1.2, which is consistent 
with X-ray temperature profiles of nearby clusters ( Markevitch 
et al.|p^98] ), and fix the normalisation of p by requiring a 
baryon fraction of 0.12 at the scale radius, which we set to 
rs = 300/z~^kpc. Note that the lengths here and below are given 
in comoving units. A reduced Hubble parameter of /z = 0.7 is 
used and c is chosen to be 1.9 x 10^ Mq / {h~^kpc)^ . These 
choices for and c correspond to a massive galaxy cluster. To 
test the reconstruction methods we put this analytically modelled 
cluster at a redshift of z = 0.3 and produce synthetic X-ray, ther- 
mal SZ and lensing observations. 



3.2. Synthetic observations 

Eor the X-ray observatio ns we use the table of the cooling func- 
tion discussed in section 2.1 and a 128 x 128 x 128 grid with 
1.5/z~^Mpc side length to project the gas distribution and get 
a map of the X-ray surface brightness in a 0.25-7.0 keV band. 
Then, except for reconstructions we specifically characterise as 
done without observational noise, we add photon noise corre- 
sponding to 10^ observed source photons to these maps using 
the same method as in Paper I. 

The thermal SZ maps are generated like in Paper I by project- 
ing the product of cluster gas density and temperature along the 
line-of-sight onto a 128 x 128 grid with 1.5/z~^Mpc side length. 
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Then the result is appropriately rescaled and, unless stated other- 
wise, noise corresponding to future ALMA Band 3 observations 
is added. In Band 3 (84-116 GHz) and in its compact configura- 
tion, ALMA will be able to achieve a temperature sensitivity of 
50/iK at a spatial resolution of ~ 3 arcsec in about four hours of 
observation (see [Butler & Wootten|1999| ). 

To produce maps of the lensing potential we project the mass 
inside a cube of 6/z~^Mpc side length which is centred on the 
cluster along the line-of- sight and calculate the convergence. 
The grid we use for this purpose is chosen such that each pixel 
corresponds to roughly 1/3 square arcminute on the sky, so that 
it contains about 10 galaxies, if we assume an average density 
of background galaxies useable for a weak lensing analysis of 
Ug = 30/arcmin^. For instance for a cluster at redshift z = 0.3 a 
44 X 44 convergence map covers the projection of the 6/z~^Mpc 
cube on the sky. 

For lensing reconstructions with observational noise, nor- 
mally distributed noise with variance 

is added to each pixel of the convergence map. Here cr^ is the 
variance expected for a weak lensing reconstruction of the con- 
vergence for a density Ug of background galaxies with an intrin- 
sic ellipticity dispersion CTg, and for an angular pixel size a of 
the convergence map (see |van Waerbeke 2000). It is assumed 
that the galaxy ellipticities are smoothed with a Gaussian of an- 
gular standard deviation G before the reconstruction. We choose 
G = a and G^ = 0.3. 

Then the convergence map is used to calculate the lens- 
ing potential in Fourier space. A source redshift of 1.5 is as- 
sumed. To reduce errors introduced by the implicit assumption 
of a periodic convergence field of such Fourier methods (see e.g. 
puchwein et al. 2005), we first zero-pad the convergence map 
to 2048 X 2048 pixels. Once the map of the lensing potential is 
calculated, we crop it to its original size before using it for re- 
constructions. 

3.3. Cumulative mass profiles 

We applied these methods to produce synthetic X-ray, SZ and 
lensing observation for the analytic halo described above. We 
then generate three-dimensional reconstructions of the cluster 
gas based on X-ray and SZ data and reconstructions of the grav- 
itational potential based on lensing data using the methods de- 
tailed in section [2] and in Paper I. 

Figure [T] shows the cumulative mass profiles M<^ xsz(^), 
obtained from the XSZ reconstructions, and M<^ iensing(^), ob- 
tained from the lensing reconstructions, and compares them to 
the original analytic profile. The profiles are shown for recon- 
structions based on data without observational noise and for 
reconstructions based on noisy data and for inclination angles 
/ = 30° and / = 70° between the symmetry axis and the line- 
of-sight. The inclination angles were assumed to be known for 
the reconstructions. See Paper I for methods to determine them 
from the observations. For the gas reconstructions ratios / of 
V/?/p to — V0 of / = 1.0 and / = 0.8 were used. The lensing 
reconstructions were done using both flat priors and spherically 
symmetrised priors. 

The XSZ and the lensing mass profiles agree very well for 
a halo in hydrostatic equilibrium (/ = 1.0) and when using data 
without noise and an inclination / = 70° (see upper right panel 
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Fig. 2. Cumulative mass profiles M<^(r) of an ellipsoidal an- 
alytic halo and its reconstructions from lensing maps without 
noise. Inclination angles of / = 30° and / = 70° were used for 
the synthetic observations and assumed to be known for the re- 
constructions. Lensing reconstructions are shown for a flat and 
for the spherically symmetrised prior. 



of Fig. [T]). They also excellently match the original analytic pro- 
file. The only significant difference between the profiles is that 
the lensing mass is too small very close to the cluster centre. 
However, this is completely expected because the lensing ob- 
servations lack the resolution required to accurately resolve this 
region. When perturbing the hydrostatic equilibrium by 20%, 
in other words when assuming / = 0.8, the XSZ reconstructed 
mass profile is essentially 0.8 times the original analytic pro- 
file as theoretically expected. For such a halo one can easily see 
a significant (~ 20%) difference between the lensing and XSZ 
mass profiles, which directly reflects the deviation from hydro- 
static equilibrium. Also when adding noise to the synthetic ob- 
servations (see lower right panel) such a deviation from hydro- 
static equilibrium can be faithfully reproduced. For smaller in- 
clination angles of / = 30° (left panels) the accuracy of the re- 
constructions is somewhat lower and one can also see significant 
differences between the lensing reconstructions based on a flat 
prior and a spherically symmetrised prior. The latter reproduce 
the original profiles much better. Thus, for such small inclina- 
tions one can detect deviations from hydrostatic equilibrium by 
comparing lensing reconstructions based on a spherically sym- 
metrised prior to XSZ reconstructions. Also note that for a ran- 
domly oriented cluster sample only about 13% of the clusters 
have inclination angles smaller than 30°. Paper I contains a more 
detailed discussion of the dependence of the reconstruction ac- 
curacy on the inclination angle. 

Above we used spherically symmetrised priors in the lensing 
reconstruction of spherically symmetric halos. It is reassuring, 
but not really surprising that this works well. We thus need to 
check whether or not a spherically symmetrised prior also im- 
proves the lensing reconstruction quality of elliptical halos for 
small inclination angles. In Figure |2] we show lensing recon- 
structions of the cumulative mass profile of an elliptic analytic 
halo with an NFW density profile but isodensity surfaces that 
are prolate spheroids with a major to minor axis ratio of 2 to 1. 
The lensing reconstructions with a spherically symmetrised prior 
reproduce the original analytic profile well, both for small and 
for large inclination angles. On the other hand when using a flat 
prior we again obtain too small lensing masses for small inclina- 
tion angles. It is thus favourable to use a spherically symmetrised 
prior for the iterative deprojection of the lensing potential. 
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Fig. 1. Cumulative mass profiles M<^(r) of an analytic halo and its reconstructions from X-ray and SZ maps with and without 
observational noise as well as from lensing maps with and without noise. The upper panels show the results obtained from maps 
without noise, while the lower panels show the profiles found from noisy maps. For the reconstructions shown in the left panels an 
inclination angle of 30° was used and assumed to be known, while the right panels show the corresponding results for an inclination 
angle of 70° . Lensing reconstructions are shown for a flat and for the spherically symmetrised prior. The XSZ reconstructions were 
done for halos with ratios / of Vp/p to —Vij) of 1.0 and 0.8. For comparison we also show the original analytic cumulative mass 
multiplied by 0.8. 



3.4. Gravitational potential reconstructions 

In Figure [3] we show three-dimensional reconstructions of the 
gravitational potential of the analytic halo from X-ray and SZ 
data and from lensing data as well as the original analytic grav- 
itational potential described by Eq. (|7|. Reconstructions that are 
based on idealised observations without noise and on more re- 
alistic noisy observations are shown. The XSZ potential recon- 
structions were obtained from the X-ray, SZ cluster gas recon- 
structions by assuming hydrostatic equilibrium and by using the 
minimisation method described in section |2.2| The lensing re- 
constructions were obtained directly by deprojecting the lensing 
potential. The XSZ reconstructions reproduce the inner region of 
the cluster well, while the lensing reconstructions lack the res- 
olution to accurately resolve this innermost part. Between dis- 
tances r from the cluster centre of 150/z~^kpc and 450/z~^kpc 
both reconstruction methods yield very good results. Farther 
outside the lensing reconstruction is still accurate, while the 
XSZ reconstruction becomes more and more unrealistic. This 
is partly due to different noise properties. But in the example 
shown in Figure [3] it is also due to the smaller box size used 
for the XSZ reconstructions and reconstruction artefacts that de- 
velop close to the perimeter of this box. The weighting func- 
tion w(r), which was introduced to prevent non-local effects of 
these artefacts, was chosen to decrease from unity to zero be- 
tween r = 450/z~^kpc and r = 600/z~^kpc in these XSZ recon- 
structions. Thus they become unrealistic farther outside. 



4. Probing the dynamical state of numerically 
simulated clusters 

4.1. The simulations 

In the previous section, novel methods to probe hydrostatic equi- 
librium in clusters of galaxies were tested with synthetic obser- 
vations of analytically modelled clusters. For a more realistic 
test we now apply these methods to a sample of four numeri- 
cally simulated galaxy clusters. The same sample was also used 
in Paper I. The simulations were carried o ut by Klaus Dolag 
with the GADGET-2 code (Springel 2005), a new version of 



the parallel TreeSPH simulation code GADGET ( [Springel et al. 
2001 ). The cluster regions were extracted from a dissipation-less 
(dark matter only) simulation with a box size of 419h~^ Mpc 
of a flat ACDM model with Q^rn = 0.3, h = 0.7, dg = 0.9 (see 
Yos hida et al.||200 1 ). They were re-simulated with higher res- 
olution using t he "Zoomed Initial Conditions" (ZIC) technique 
( Tormen|1997 i Gas was introduced into the high-resolution re- 
gion by splitting each parent particle into a gas and a dark mat- 
ter particle, which were then displaced by half the mean inter- 
particle distance, such that the centre-of-mass and the momen- 
tum were conserved. The mass ratio of gas to dark matter par- 
ticles was set to obtain Q.^ = 0.04. The final mass resolution 
was rriDM = 1.13 x 10^ h~^MQ and mgas = 1.1 X 10^ h-^MQ for 
dark-matter and gas particles within the high-resolution region, 
respectively. The simulations we use follow the dynamics of the 
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Fig. 3. Gravitational potential of an analytic halo and its reconstructions from synthetic X-ray and SZ observations and from syn- 
thetic lensing observations each with and without observational noise. The XSZ reconstructions work well even close to the cluster 
centre, where lensing observations lack the resolution to accurately resolve the central peak. However, farther outside the lensing 
reconstructions perform better, due to the different noise properties, but in the example shown here also because of the smaller box 
size used for the XSZ reconstructions. 



dark matter and the adiabatic evolution of the cluster gas, but 
they ignore radiative cooling. They are described in more detail 
in |Puchweinetar| ( [2005] ) and |Dolag et"aL] ([20Q5). 



4.2. Cumulative mass profiles and hydrostatic equilibrium in 
simulated clusters 

We produced synthetic X-ray, thermal SZ and lensing observa- 
tions of these four simulated clusters for 28 simulation snapshots 
between redshifts 0.58 and 0.1 and threeUnes -of- sight using es- 
sentially the same methods as in section 3.2 for the analytic halo. 
The only difference is that we did not use a three-dimensional 
grid for projections along the line-of-sight. For the X-ray and 
SZ maps the X-ray luminosities and integrated Compton y pa- 
rameters of the gas particles are projected directly onto a two- 
dimensional 128 X 128 grid by using the particles' projected 
SPH smoothing kernel. The convergence of the simulated clus- 
ters is found in a similar way by projecting the masses of both 
gas and dark matter particles onto a two-dimensional grid, whose 
dimensions are again chosen such that one pixel corresponds to 
roughly 1/3 square arcminute on the sky. O bser vational noise is 
added in exactly the same way as in section [3]2] 

Based on these synthetic observations we perform three- 
dimensional XSZ reconstructions of the cluster gas distribution 
and lensing reconstructions of the gravitational potential. The in- 
clination angle is assumed to be known for the reconstructions. 
Paper I discusses how it can be determined from data and how a 
symmetry axis is chosen for the simulated clusters. Spherically 
symmetrised priors are used for the lensing reconstructions. The 
reconstructions are then used to probe hydrostatic equilibrium 
by calculating and comparing their cumulative mass profiles 

M<r, XSZ (r) and M<,, lensing {r) . 

In Figure |4] we show these profiles for two clusters that did 
not experience a major merger recently. For comparison we also 
show the original simulated mass profile and the profile that 
would be expected from the original simulated gas distribution 
by assuming hydrostatic equilibrium. The latter is calculated like 
^<r,xsz(^), however directly from the simulated gas distribution 
rather than the reconstructed one. We again use 128 rings which 
are equally spaced in the polar angle 6 to numerically evaluate 
the surface integral in Eq. Q, but as the simulated gas distribu- 



tion is not perfectly axisymmetric we use 128 sampling points 
equally spaced in the longitude angle for each of these rings. 
The gas density p and the pressure gradient Vp are calculated at 
each sampling point using the SPH formalism, i.e. by summing 
up the contributions from all nearby particles using their SPH 
smoothing kernels and the gradients thereof. For these relaxed 
clusters the XSZ reconstructed profiles and the lensing recon- 
structed profiles agree well with each other, with the original 
mass profile and the profile obtained from the original gas dis- 
tribution. This shows that for such relaxed clusters this method 
allows accurate and consistent lensing and XSZ mass estimates. 
The results also confirm that these clusters are close to hydro- 
static equilibrium. 

It is reassuring that this novel method to probe hydrostatic 
equilibrium works well for clusters that do not have a record of 
recent mergers. However clusters that do experience such vio- 
lent events may be even more interesting to study. In Figure |5] 
we show a cluster at four different times during a merger. For 
each of these snapshots we show reconstructions of the cumula- 
tive mass profile from synthetic X-ray, SZ and lensing observa- 
tions, as well as the original mass profile and the profile obtained 
from the simulated gas distribution. Again observational noise 
was added to the synthetic maps used for the reconstruction. We 
also show X-ray maps of the cluster for each of the four snap- 
shots. These are however idealised noise-free versions and just 
meant to illustrate what is going on in the cluster. To facilitate 
following the merger we also show the approximate trajectory 
of the relevant inf ailing subhalo in the X-ray maps. 

In the first snapshot (upper left panel) the main cluster halo is 
still close to hydrostatic equilibrium. The lensing and XSZ mass 
estimates still agree well for radii r smaller than the distance to 
the inf ailing subhalo. In the second snapshot (upper right panel), 
after the subhalo has passed the main halo, shocked gas causes a 
too large XSZ mass estimate from roughly the subhalo distance 
outwards. The mass profile obtained directly from the simulated 
gas distribution shows the same behaviour and thus confirms that 
this is not an artefact of the reconstruction but a real, signifi- 
cant deviation from hydrostatic equilibrium, which is recovered 
by the reconstruction or in this example even somewhat over- 
estimated. The lensing reconstruction still reproduces the origi- 
nal simulated mass profile well. Thus by comparing lensing and 
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Fig. 4. Cumulative mass profiles M<^(r) of relaxed simulated clusters gl at redshift z = 0.25 and g51 at redshift z = 0.3. Profiles 
of the original simulated mass distribution, of the lensing and of the XSZ reconstructions are shown, as well as the profile obtained 
directly from the simulated gas distribution by assuming hydrostatic equilibrium. The lensing and XSZ reconstructions are based 
on synthetic observation that contain observational noise. For such relaxed clusters both the lensing and the XSZ reconstructions 
agree very well with the original mass profile. 



XSZ cumulative mass profiles one can directly see the devia- 
tions from hydrostatic equilibrium. The third snapshot (lower 
left panel) shows that when the bow shock moves outward one 
can also obtain too low cluster masses by assuming hydrostatic 
equilibrium during a merger. Again the effect can be seen in both 
the mass profiles obtained directly from the simulated gas distri- 
bution and obtained from the three-dimensional XSZ gas recon- 
struction. In the fourth snapshot (lower right panel) hydrostatic 
equilibrium is already largely restored, even if one can still see 
the pronounced bow shock in the X-ray map. 

These simulations show that deviations from hydrostatic 
equilibrium during mergers can be faithfully recovered by the 
cluster reconstructions methods introduced above. 

4.3. Accuracy and reliability of cumulative mass profile 
reconstructions 

To determine the typical scatter in cumulative mass profile re- 
constructions and quantify the significance of detections of devi- 
ations from hydrostatic equilibrium we repeated the reconstruc- 
tion of the merging simulated cluster shown in the upper right 
panel of Fig. [5] with different noise realisations and for different 
lines-of-sight. 

For the left panel of Fig. [6] we used the same line-of-sight as 
in Fig. [5] but 50 different noise realisations for the synthetic X- 
ray, thermal SZ and lensing observations. The noise realisations 
were obtained using different seeds for the random number gen- 
erator employed for adding noise to the synthetic observations. 
The mean XSZ and lensing reconstructed profiles and the l-cr 
errors are shown as well as the profile of the original simulated 
mass distribution and the profile obtained directly from the sim- 
ulated gas distribution by assuming hydrostatic equilibrium. The 
deviations from hydrostatic equilibrium are reliably detected. As 
expected for a cluster that contains substructure that violates ax- 
ial symmetry there are also some systematic deviations such that 
the mean profiles are not centred exactly on the simulated pro- 
files. 

For the right panel we started with a sample of synthetic lens- 
ing. X-ray and SZ observations along 50 different randomly ori- 



ented lines-of-sight. All contain realistic observational noise. It 
turned out that for projections for which the merging subhalo 
responsible for perturbing hydrostatic equilibrium is almost di- 
rectly in front of or behind the main halo detecting deviations 
from hydrostatic equilibrium is less reliable. This is not surpris- 
ing as the signal from the region where hydrostatic equilibrium is 
strongly perturbed is superimposed with a larger signal from the 
main halo, so that the contributions to such projections are diffi- 
cult to separate. For the right panel of Fig.|6]we thus decided to 
reject all 16 lines-of-sight for which the projected distance of the 
relevant subhalo from the main halo centre is less than 200 
kpc, as well as one line-of-sight which happened to be inclined 
by only 2° with respect to the cluster's symmetry axis which 
is to small for a faithful reconstruction. The mean and the l-d 
errors of the reconstructions that were based on the 33 remain- 
ing lines-of-sight are shown. Again deviations from hydrostatic 
equilibrium can be reliably detected. 

As discussed in Paper I reconstruction artefacts can appear 
close to the perimeter of the box used for the reconstructions. 
As we can see in the right panel of Fig. [6] they can dominate the 
XSZ reconstructed cumulative mass profiles errors from roughly 
r = 600h~^ kpc outwards for some lines-of-sight, when using a 
1.5h~^ Mpc sidelength box for the XSZ reconstruction. Thus 
when the quality of the observations allows studying a larger 
region one should also use an appropriately larger box for the 
XSZ reconstruction to avoid this problem. 



5. Summary and discussion 

We proposed a novel method to obtain three-dimensional recon- 
structions of a galaxy cluster's gravitational potential and cumu- 
lative mass profile from X-ray and thermal SZ observations un- 
der the assumption of hydrostatic equilibrium and independently 
dropping this assumption from lensing data. If only X-ray and 
thermal SZ data is available accurate reconstructions of relaxed 
clusters can be obtained. If, however, lensing data is available as 
well, hydrostatic equilibrium can be probed, also in dynamically 
active clusters, by comparing these independent reconstructions. 
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Fig. 5. Cumulative mass profiles M<^(r) and X-ray surface brightness maps of simulated cluster gSlat four different redshifts during 
a merger. The approximate trajectory of the inf ailing subhalo is illustrated in the X-ray maps. Profiles of the original simulated mass 
distribution, of the lensing and of the XSZ reconstructions are shown, as well as the profile obtained directly from the simulated 
gas distribution by assuming hydrostatic equilibrium. The lensing and XSZ reconstructions are based on synthetic observation that 
contain observational noise. The X-ray maps shown above are however idealised noise-free versions and were rotated such as to all 
have the same orientation in space. 
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Fig. 6. Mean XSZ and lensing reconstructed cumulative mass profiles M<^(r) and their 1-a errors of merging simulated cluster 
g51 obtained for different noise realisations (left panel) and different lines-of-sight (right panel). Profiles of the original simulated 
mass distribution and the profiles obtained directly from the simulated gas distribution by assuming hydrostatic equilibrium are 
shown for reference. X-ray, SZ and lensing observations along one line-of-sight but with 50 different noise realisations were used 
for the reconstructions whose mean and 1-cr errors are shown in the left panel. For the right panel we started with a sample of noisy 
synthetic observations along 50 different randomly oriented lines-of-sight. However we rejected projections for which the merging 
subhalo is almost directly behind or in front of the main halo (projected distance < 200 kpc) as well as one line-of-sight with 
an inclination of only 2° with respect to the clusters symmetry axis which is to small for a faithful reconstruction. The mean and the 
l-cr errors of the profiles reconstructed from the observations along the remaining 33 lines-of-sight are shown. For both the different 
noise realisations and the different lines-of-sight deviations from hydrostatic equilibrium can be reliably detected. 



The three-dimensional reconstructions are based on iterative 
Richardson-Lucy deconvolution and assume only axial symme- 
try of the cluster halo with respect to an arbitrarily inclined axis. 
The X-ray and thermal SZ data are used to first reconstruct the 
three-dimensional cluster gas density and temperature distribu- 
tion. No equilibrium assumption except local thermal equilib- 
rium is needed for that. Then the gravitational potential and the 
cumulative mass profile can be obtained from these reconstruc- 
tions under the assumption of hydrostatic equilibrium. For the 
lensing reconstructions we deproject the lensing potential ob- 
tained by a weak lensing or a combined weak and strong lensing 
analysis. This yields the three-dimensional gravitational poten- 
tial, from which we can get independent cumulative mass pro- 
files by exploiting Gauss's law. The X-ray and thermal SZ anal- 
ysis (abbreviated by XSZ throughout this work) and the lensing 
analysis are then compared in order to probe hydrostatic equi- 
librium and to test the accuracy of mass estimates based on the 
assumption of hydrostatic equilibrium. 

These methods were tested with synthetic X-ray, thermal SZ 
and lensing observations of analytically modelled and numeri- 
cally simulated galaxy clusters. Except where specifically noted 
realistic observational noise was added to the synthetic observa- 
tions. 

For analytically modelled clusters in hydrostatic equilibrium 
we found: 

- Consistent and accurate lensing and X-ray, SZ based cumu- 
lative mass profiles lensing(^) and M<^ xsz(^) can be ob- 
tained. 

- The accuracy somewhat decreases for very small inclination 
angle between the line-of-sight and the cluster's symmetry 
axis. 

- Higher accuracy for the iterative deprojection of the lens- 
ing potential for small inclination angles are achieved with 
spherically symmetrised priors. 



- Faithful three-dimensional reconstructions of the gravita- 
tional potential can be obtained from both lensing observa- 
tions and from an XSZ analysis. 

For analytically modelled clusters that are not in hydrostatic 
equilibrium, we showed that the deviations from equilibrium can 
be effectively probed by a comparison of lensing and XSZ recon- 
structions even when realistic observational noise is present. 

From reconstructions based on synthetic observations of a 
sample of numerically simulated galaxy clusters we conclude: 

- For clusters that did not experience recent mergers consistent 
and accurate lensing and XSZ cumulative mass profiles are 
found. 

- Although these clusters are not perfectly axisymmetric and 
noise is added to the synthetic data, the accuracy of recon- 
structed cumulative mass profiles is typically better than 10 
to 15% for both the X-ray, SZ and the lensing reconstruc- 
tions. 

- On the other hand in clusters in the process of merging devia- 
tions from hydrostatic equilibrium can be accurately probed, 
except for cases where the relevant merging subhalo appears 
directly in front of or behind the main halo's centre. 
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